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1. Introduction 

First principle predictions for the QCD phase diagram would be of great value as a benchmark 
for the international experimental heavy ion program. In principle, we know precisely how to 
proceed: start from the grand canonical partition function, 

J dAdet 2 (D + HYo + m)e' SYM , (1.1) 

and study the baryon number as a function of the chemical potential, [i, and temperature for quark 
masses m (for notational simplicity we consider the two flavor theory in these proceedings). As the 
phase transitions occur in the non perturbative domain it is natural to turn to lattice QCD. While 
we know exactly how to include the chemical potential on the lattice [jl|, g], in practice, our studies 
are limited by the fact that the fermion determinant at nonzero values of the chemical potential 
becomes complex 

det 2 {D + HY + m) = |det(D + iiy {) + m)\ 2 e 2ie . (1.2) 

This sign problem invalidates the standard Monte Carlo method which is at the hart of lattice QCD. 
The sign problem is not only hard because the average of the phase factor is exponentially small 
[gl, it is also challenging because much of our intuition for statistical systems leads to wrong 
conclusions at non zero ii. Most prominently, a probabilistic argument leads one to conclude that 
the chiral condensate is continuous as a function of the quark mass at m = for any non zero 
value of ii [f|]. Rather, as the solution [j|] of this Silver Blaze problem shows, it is imperative 
that we consider distribution functions that take complex values. Only through extreme complex 
oscillations of the eigenvalue density of the Dirac operator is it possible to understand how the 
discontinuity of the chiral condensate remains non zero in the presence of the chemical potential. 
Here we will show [Q] that the extreme complex oscillations take place also in the chiral condensate 
and the baryon number and that they are essential to get the correct physical results. 

For simplicity we will here focus on the baryon number, which in a single gauge field config- 
uration is given by 

n = — log det(D + jljo + m). (1.3) 

Our goal is to determine the distribution (8(n — n')) of the baryon number over the gauge fields. We 
will show that this distribution takes complex values and that the extreme oscillations are essential 
to obtain the correct average baryon number 

(n) = [ dn n (S(n-n)}. (1.4) 



This direct insight into the sign problem also allow us to address how complex Langevin works in 
this case. 

In order to understand better the distribution (8(n — n')) of the baryon number operator over 
the gauge fields, let us first understand its first two moments. For the first moment there is no source 
of confusion: The average quark number is the first moment 

b^ Zi «=<">- (15) 
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However, the second derivative with respect to /I 

1 1 d 2 



4Z^^ 1+1 = ( " 2) + K^)' 



(1.6) 



is not the second moment of the distribution (8 (n — n')). Rather, the average of the square of n can 
be written as 



(n 2 ) 



1 d d 
Z d\i u djld 



Zi+i 



(1.7) 



where we distinguished the chemical potentials for the two flavors and differentiated with respect 
to each one before setting them equal. 

When we express the traces in terms of the eigenvalues, Zk, of Yq(D + m) 



1 1 



2Zi + i d\l 
1 1 d 2 



4+1 



4Zi + i d\l 
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Y % + M / 
1 1 



(1.8) 
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it becomes obvious that ( |L~6| ) is not the average of a square and in particular it is not the second 
moment of the distribution of n over the gauge fields. The distribution (8(n — n')} is nevertheless 
of great interest since it gives direct insights in the sign problem. 

As a final point before we turn to the results, note that the quark number takes complex values 



Tr- 



7o 



n(-n). 



Hence, the distribution (8(n — n')) is in the complex n plane 

P n (x,y) = (8 (x - Re[n]) 8 (y - Im[n])) , 



(1.9) 



(1.10) 



and the average baryon number is given by the integral of (x + iy) weighted by the distribution 

P n (x,y)dxdy. 



2. The distribution of n from Chiral Perturbation Theory 

Despite the fact that pions have zero baryon charge the distribution of the baryon number 
over the gauge fields is non trivial when computed within Chiral Perturbation Theory. Certainly in 
Chiral Perturbation Theory we have that 
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Figure 1: Scatter plot of the spectrum of a Random Matrix Dirac operator for jj. — m K /2; left: the eigenval- 
ues of D + jiiyo + tn right: the eigenvalues of i(yb(D + m) + jit). In both cases the support of the spectrum 
has reached the origin indicated by the red point. Beyond this point, ie. for /J. > m x /2 the distribution of 
the chiral condensate and the baryon has power law tails. (A similar phenomenon is expected to happen for 
lattice QCD with Wilson fermions in the Aoki phase [||.) 



but the average of the square of n is non zero 

#0, (2.2) 



d 2 



d\L\d\i 2 
since the 1-loop free energy is 



GbOii,*) = I ^cosh(^,). (2.3) 

n=\ 



So (8(n — n')} must necessarily be non trivial in Chiral Perturbation Theory. 

In order to compute the full distribution P n (x,y) it is necessary to evaluate all moments 

(Re[n] k lm[n] j ) (2.4) 

in Chiral Perturbation Theory. The details are given in ^ and involve an interesting combinatorial 
use of the replica trick [^]. For the computation it is essential to specify whether the chemical 
potential is larger or smaller than m w /2, since in the replicated generating functions additional 
condensates appear at this scale. Also from the perspective of the eigenvalues of the Dirac operator 
it is clear that the case }X < m n /2 must be very different from the one with }X > m n j2, see figure []]. 



We first discuss the distribution of n for \x < m n /2. To one-loop order in Chiral Perturbation 
Theory the distribution factorizes [^] 

P n (x,y) = P Re \ n ] (x)P ImW (y), (2.5) 
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Figure 2: The distribution of the baryon number over the gauge fields in the grand canonical ensemble. 
For ji < m n /2 the distribution factorizes into the distribution of the real part (left figure) and the imaginary 
part (right figure). The distribution of the imaginary part takes complex values, shown is the real part. Note 
the difference in the scales on the vertical axis. The amplitude of the distribution of the imaginary part of n 
grows exponentially with the volume. 



where the two factors take simple Gaussian forms 

1 

*Re[n]W = 
P\m[n](y) = 



i 

< 



-(x-Vi) 2 l(X B ull +x' ud ) 



-(y-iv,) 2 /{xL-x B ud ) 



(2.6) 
(2.7) 



Note in particular that P\ m ^ (y) takes complex values. It is quite natural that the sign problem 
manifest it self in the distribution of the imaginary part of n since 



n = -^logdet(D + ^7o + m) = -^-log|det(D + i U7o + m)| 

In the above expression for the distribution of n we have made use of the notation 

d 



(2.8) 



d}X\ 



(2.9) 



Xud 



Xud 



d\l\d\l 2 

d 2 
d\l\d\i 2 



AGoGlii,^ 
AG (-jUi,/i2) 



where the free energy difference is (note that % l ud + Xud > an d xid ~ Xud > 0) 



.mlT 2 



cosh( — - — n) - 1 



(2.10) 



Since all of the above quantities are extensive the amplitude of P Im [„] (y) grows exponentially with 
the volume, the width grows like y/V, while the period of the oscillations are of order V°. For a 
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plot of the distribution see figure ^. The extreme oscillations of P^y (y) are essential in order to 
obtain zero expectation value of the quark number in Chiral Perturbation Theory 

(n) = J dxdy (x + iy)P n (x,y) (2.11) 

= Jdx xP Re[n] (x) + ijdy yP Im[n] (y) 
= V/ + iiVj = 

The detailed cancellation between the contribution from the real part and the imaginary part is only 
possible if the phase of the fermion the determinant is accounted for properly. Similar cancellations 
also take part for the higher moments of the baryon number as well as for the moments of the chiral 
condensate, see [||] for details. 



3. Complex Langevin 

One use of the results for the distribution of the baryon number is to illustrate how the complex 
Langevin method can deal with sign problems in simple models. Clearly the distribution of the 
imaginary part of the baryon number over the gauge fields is the challenging part. Therefore, let us 
ask if complex Langevin is able to do the simple one dimensional integral 



J dyy Pjm[„](y), (3.1) 



that is, to measure the contribution to the average baryon number from the imaginary part of n. 
To this end we define the complex Langevin action for y = lm[n] as 

S = -log[P lm[n] (y)} = -{iy + Vl f/{ x l d - X f id ). (3.2) 

The next step is to complexify \m[n] asy = a + ib and write down the flow equations for a and b 

a n+ \ = a n — e— — — g- + Vsrj n (3.3) 

%ud %ud 

2(b n - v 7 ) 
Pii+l - o n - e— -j _ B . 

&ud A u d 

Note that the flow equations decouple. The equation for a is that of a Gaussian for which complex 
Langevin works perfectly. That of b simply shifts y by V/ in the imaginary direction. Since there 
is no noise in the imaginary direction, the complex Langevin method effectively shifts the contour 
of the ^-integral by a term of order V in the imaginary direction. After the shift, a simple integral 
over a Gaussian without oscillations is left and the complex Langevin method has no problem in 
evaluating this. Clearly the shift of the contour is the only reasonable thing to do in this case, the 
strength of the complex Langevin method is that it can make this shift automatically. A similar 
example was worked out in ^ and JT0|]. 

For jU > m n /2 the chemical potential enters the spectral support of /o(D + ra) and the distri- 
bution of the baryon number develops power law tails [Hj. Nevertheless, complex Langevin is also 



able to deal with the sign problem for one dimensional QCD [11] in this region. 
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4. Summary 

The interplay between lattice QCD and analytical studies of QCD is essential to understand 
QCD at nonzero chemical potential. Due to the sign problem, the standard methods of lattice 
QCD only have a limited range of applicability. In order to study dense strongly interacting matter 
from first principles new numerical methods must be invented and put to use. To understand how 
such methods can be designed it is essential to understand how the sign problem affects physical 
observables such as the baryon number and the chiral condensate. 

Here we have derived the distribution of the baryon number over the gauge fields from Chiral 
Perturbation Theory. We have shown that the distribution takes complex values and is strongly 
oscillating. These oscillations were shown to be central to the detailed cancellations which take 
place when forming the average baryon number. The distributions also give detailed information 



on the overlap problem as will be discussed in []12|]. Here we have used the distribution of the 



baryon number to show how the complex Langevin method can deal with sign problems. An 
important point to take away from this is that the complex Langevin method works equally well 
independent of the volume and hence independently of the strength of the sign problem. Similarly, 



in the well known cases [13] where the complex Langevin method fails it does so independently of 
the volume. 

It is also possible to compute the distribution of the baryon number over the phase of the 
fermion determinant within Chiral Perturbation theory Also in this case the complex and 

oscillating nature of the distribution is essential in order to obtain the correct physics at nonzero 
chemical potential. That calculation also directly demonstrates that all phases of the fermion deter- 
minant are important. 
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